Stationary and moving breathers in a simplified 
model of curved alpha— helix proteins 



JFR Archilla f , Yu B Gaididei }, PL Christiansen § and J 
Cuevas f 

fNonlinear Physics Group of the University of Sevilla. Dep. Fi'sica Aplicada I, 
ETSI Informatica. Avda. Reina Mercedes s/n, 41012 Sevilla, Spain. 
JBogolyubov Institute for Theoretical Physics. 03143 Kiev, Ukraine. 
§Informatics and Mathematical Modelling. 

The Technical University of Denmark, DK-2800 Lyngby, Denmark. 
E-mail: archilla9us.es 



Abstract. The existence, stability and movability of breathers in a model for 
alpha-helix proteins is studied. This model basically consists a chain of dipolc 
moments parallel to it. The existence of localized linear modes brings about 
that the system has a characteristic frequency, which depends on the curvature 
of the chain. Hard breathers are stable, while soft ones experiment subharmonic 
instabilities that preserve, however the localization. Moving breathers can travel 
across the bending point for small curvature and are reflected when it is increased. 
No trapping of breathers takes place. 
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1. Introduction 

In the last years quite a great deal of research is being done on localization due to the 
interplay of nonlinearity and geometry, either with the Nonlinear Schrodinger Equation 
§ |], FPU models {§|, §, Klein-Gordon models §, § or in DNA models g §. The 
objective is to understand the role of the bending points in biomolecules: whether 
localized excitations can travel across them or not and whether there are points where 
energy is stored and play a biological function. The description of the biological 
systems is given by variables, that represent internal or external degrees of freedom, 
which oscillate with time and are coupled by different potentials. A change of geometry 
can be felt by the system by different physical mechanisms, modelled correspondingly: 
interaction between nearest and next neighbours [Q, ; potentials that depend on the 
angles or long-range interaction due to the dipole-dipole coupling [§, @, [|. 

Simplifying, the effect of a curved chain can be described easily. For stationary 
excitations the zone where the chain of oscillators is bent is inhomogeneous, bringing 
about the existence of localized linear modes which compete strongly with the 
nonlinear localized modes. For high coupling the linear localization predominates, 
for low coupling the nonlinear localization is the important one. The transition from 
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one regime to the other can be continuous or discontinuous as the parameters are 
changed or the curvature increases Moving excitations, have to travel across this 
inhomogeneous zone. The excitation within the bending zone needs a different energy, 
which can be larger than in the straight system, therefore acting as a barrier. The 
excitation can be reflected, transmitted or trapped pC| |. 

The alpha-helix protein is another molecule for which is interesting to investigate 
the role of the bending in the existence, shape, properties and transport of localized 
excitations. The peptide groups have a dipole moment parallel to the chain and the 
Amide-I excitations interact among them through acoustic phonons. The dipole-dipolc 
interaction to be described in detail below makes the system sensitive to the shape 
of the molecule. Apart from the interest of the biophysical problem in itself, it has 
a relevant feature from a more theoretical point of view: it is only necessary to take 
into account the nearest neighbours in order that the system can feel the bending of 
the chain. To our knowledge the effect of the curvature in such a system with dipoles 
parallel to the chain have not yet been described. 

There are different excitations that can be considered, as solitons, envelope 
solitons and discrete breathers. Here we focus our research on the latter. As it is well 
known, discrete breathers are very localized oscillations that appear as a consequence 
of nonlinearity and discreteness |ll], |l2, 13, 14J. They are, therefore, specially relevant 



in biomolccules when considering excitations that involve only a few units, far from the 
continuous limits. Although a relatively new field, stationary breathers are now well 
understood and research is presently focusing on the possible physical and biological 
consequences of their existence and in future technological applications, specially with 
Josephson-junctions |l6). Moving breathers is not such a mature field: there are 
techniques to obtain them but many questions remain unanswered [ p"7[ |l8| ]. Perhaps 
the most important one is the reason why they exist in some systems and not in others. 
New developments are given in Q • 

This study intends to complete previous works || [j], [nj on the properties of 
localization and transmission of energy in bent chains, described by Klein-Gordon 
models. The models in the references cited are inspired in DNA, with three main 
properties: 

(i) There is long-range interaction between dipole moments. 

(ii) The dipoles moments are perpendicular to the chain and to the plane of curvature. 

(iii) The effect of the curvature enters the model as a shortening of the distances 
between dipole moments. 

The simplified alpha-helix model studied here, apart from its different physical 
origin, has the following differences: 

(i) The only interactions considered are nearest-neighbour. Note that the DNA 
models cited, without long-range interaction would not feel at all the shape of 
the molecule. 

(ii) The dipole moments are parallel to chain and to the plane of curvature 

(iii) The bending of the chain appear in the model as a change in the angles between 
dipole moments, and, therefore in the interaction energy. 

In spite of these differences, both models present very similar behaviour both for 
predominant stronger attractive or repulsive interaction. This suggest that they are 
generic on bent Klein-Gordon systems. 
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However, this paper introduces a new point of view. Previous works 
misinterpreted the mathematical property of annihilation of breathers as the curvature 
is increased. What happens is that the bending forces the breather at the bending 
point to choose some frequency, which depends weakly on its energy. 

This phenomenon, could be, in principle, tested experimentally. For example, the 
mean curvature of DNA depends on the solvents concentration, and the absorbtion of 
radiation by breathers at the bending points would appear at some characteristics 
frequencies, that change with the solvents concentration. The presence of some 
harmonics of these frequencies would differentiate breathers from linear localized 
modes. Certainly, many technical problems should be expected, but we are planning 
to perform this study with some experimental groups. Apparently, the more adequate 
molecule might be RNA due to its short persistence length. 

The short section ^ deals with moving breathers. Its first objective is to check 
if in this model, breathers can move and are reflected by or transmitted through 
the bending point. The second is to confirm the trapping hypothesis introduced in 
Ref. for a DNA model with an impurity. The role of an impurity is similar 
to a bending point in bringing about localized linear modes, which compete with 
the (nonlinear) breather localization. This hypothesis consist in the nonexistence of 
trapping for moving breathers when there exists a linear localized mode with different 
tail profile, i.e., neighbouring sites in phase or with opposite phase (and, therefore, 
different frequency). There exist no mathematical proof, but it is also confirmed in 
this work in spite of the differences stated above. 



2. Model description 



We describe the position of each amino acid by r„, n being an index. The distance 
between amino acids is supposed to be a constant a, i.e. , a = |r n +i — r n |. The 

( r +1 ~ r ) 

direction of the dipole moments are given by unit vectors t„ = - and the 

| r n+l — r n 

moments themselves are given by p„ = p n t n . 

The Hamiltonian of the system, for which details are given in Appendix A, can 
be written as 

H = ^ \^il 2 n + ^o u n + *W + \^{Un+\ - U n ) 2 + /i(t„ + i - t n ) 2 U n U n+1 , (1) 

71 

with u = 1, and the on site potential is given by V(u n ) = + $(u n ) , ^f(u n ) 

being its nonlinear part. 

The corresponding dynamical equations are 

ii n + uj 2 u n + + e(2 u n - u n+ i - u„_i) 

+ /i((t„ + i - t n ) 2 u n+1 + (t„_i - t„) 2 u„_i) = , (2) 

or, in a simplified notation 

f n (u) =U n + L0 2 U n + #'(lt„) + E ^ Cn, m U m + M J ™,mU m = , (3) 

tci m 

with the obvious definition of the coupling matrices C and J. Note that C does not 
depend on the parameters while J depends only on the curvature k 

Certainly, many different shapes can be considered. Here we focus as in previous 
papers on a parabola in a two-dimensional space. The reason is that it is the 
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Figure 1. Sketch of the model for curvature K = 0.5. 

simplest form to describe a bending point and it is an approximation of any curve 
in the neighbourhood of it. Other shapes considered in the literature like the hairpin 
geometry, although interesting in themselves are not so adequate for our purposes, 
because, for example, the latter is composed by three homogeneous regions (except 
if considering long-range interaction). Therefore, the vectors r„ have components 
(xn,Un) with y n — \nx 2 n , k being the curvature. 
A sketch of the model is shown if figure [l] 

In spite of having many dimensionless variables, there are still a few parameters 
in our problem: e, which will be called the stacking coupling parameter, /i, the dipolc 
coupling parameter, the curvature k and the breather frequency, which we represent 
by LOh- The distance of lo^ from the phonon band is a measure of the degree of 
nonlinearity of the excitation. Generally speaking, the more nonlinear the excitation, 
the narrower it becomes. To explore thoroughly the parameter space is a daunting 
task but to obtain a general picture is a much more affordable one and this is what 
we intend here. 

3. Linear modes 

In bent systems there exist linear modes that are localized around the bending point. 
The first step is therefore to obtain the dependence of the linear spectrum on the 
parameters and identify the most important ones. The linear modes are the solutions 
of equation || with ^ = 0. We first obtain the dependence on the stacking parameter 
e without dipole coupling or equivalently for a straight chain, i.e., either /i = or 
(t n +i — t n ) = 0, Vrt in Equation [| Figure || shows this spectrum dependence. 

This a well known spectrum, which we include here for comparison. Let us 
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Figure 2. Dependence of the linear spectrum of the system on the stacking 
coupling parameter e for a straight chain. The insets show the shape of the linear 
modes with highest and lowest frequency. 



comment three facts: first, the spectrum is continuous; second, the mode with highest 
frequency consist of all oscillators vibrating with opposite phase with respect to the 
nearest neighbours; third, the mode with lowest frequency consist of all oscillators 
vibrating in phase. This modes will hereafter be denoted the top and the bottom 
mode, respectively. 

When the stacking coupling e = we have to choose either to represent the 
variation of the linear spectrum as a function of the parameter /i or to the curvature 
re, the other parameters fixed at non-zero values. It turns out that both possibilities 
give similar results. Figure |^ shows the dependence of the linear spectrum on re. 
There are three key features: first, there are localized modes that separate from the 
continuous spectrum; second, the top and bottom modes are localized within a radius 
of about a few units around the bending point and they become more localized as re 
increases; third, the top mode has all the oscillators in phase, which we will describe as 
being bell shaped, while the bottom mode has neighbouring oscillators with opposite 
phase (zig-zag shaped). This change of profile is due to a change of the predominant 
coupling interaction when re (or /i) is increased. Around the bending point the sign in 
front of the variable u m in the coupling terms of equation (^) changes from negative, 
i.e, attractive interaction to positive, i.e., repulsive interaction. For homogeneous 
repulsive interaction the linear spectrum would be as in figure ^ but the frequencies 
spreading downwards and the zig-zag mode at the bottom. 

The localization does not depend on the number of sites in the system. If the 
dipole coupling parameter [i is changed at constant curvature re, the spectrum is very 
similar with a slightly larger spread of the frequencies except for the top and bottom 
modes. 
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Figure 3. Dependence of the linear spectrum of the system on the curvature k 
with stacking coupling parameter e = and dipole coupling parameter fi = 0.2. 
The insets show the bell and zig-zag shapes of the linear modes with highest and 
lowest frequency, respectively. 



To conclude the description of the linear spectra we represent in figure ^ the 
dependence of the spectrum on the dipole parameter \i for a non-zero, constant value 
of the stacking parameter e = 0.05 and constant curvature k = 1. It can be seen 
that for some value of \i the same localized modes separate from the continuous 
spectrum. Note that while studying breathers these modes are going to be the ones 
competing with nonlinear excitations and be predominant for high enough coupling. 
The frequency lu^ of the breathers have to be outside the continuous band or it will 
not be possible to obtain the breathers, which in physical terms means that the linear 
localized modes will resonate with the breather frequency and the energy will spread 
along the chain. If the on-site potential is hard, the breather frequency oj-t, will be 
above the phonon band, u>h > WOj therefore competing with the bell-shaped, top, 
linear mode. In the opposite case, Wb < wo, the competing mode will be the zig-zag 
shaped, bottom one. This will to be described in detail in the next section. 

4. Hard breathers 

To obtain breathers we need to choose the value of their frequency, which for a hard 
potential (here V(u n ) — ^o u n + 1/4 ) wu l be above ujo = 1. Let us choose 
ajb = 1-2, i.e., a frequency not too far from the phonon band but outside of it. The 
properties of the breathers with stacking coupling are well known, but we include them 
here for comparison. Figure || shows the profile of the breather while the parameter e 
is changed. The path continuation using the Newton method finishes as the phonon 
band expands and collides with the breather frequency. The breather profile is a 




Figure 4. Dependence of the linear spectrum of the system on the dipole coupling 
parameter /i for constant stacking parameter e = 0.05 and curvature k = 1. The 
insets show the shape of the linear modes with highest and lowest frequency, 
respectively. 




Figure 5. Dependence of the breather profile with hard, on-site potential for 
the straight chain on the stacking parameter e. Frequency lu^, = 1.2. 



zig-zag one, which can also be easily deduced by tail analysis. 
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Figure 6. Dependence of the breather profile with hard on-site potential for 
a curved chain with curvature k = 2 on the dipole coupling parameter fi with 
constant frequency uj, = 1.2. e = 



We obtain an analogous picture if the continuation is done with respect to the 
dipole coupling parameter for constant curvature, until a bifurcation point. Now 
the breather gets the bell profile shown in figure ^| 

The fact that the breather amplitudes, i.e, {|m„(0)|} = {maxvt(|u„(i)|)}, tend 
to zero is misleading. We are in presence of an annihilation bifurcation |20| , i.e., 
the Jacobian of the dynamical equation ([3]) with respect to u = (tti, U2, ■ ■ ■ , u n ) has 
an eigenvalue that tends to zero. This eigenvalue does not correspond to another 
breather in the neighbourhood of the parameter space but to the upper (for a hard 
potential) breather band jl3| of the same breather. The physical meaning is that 
there are no localized excitation with the chosen frequency, except in the trivial case 
where the breather has zero amplitude, i.e., it is in the linear regime. We can obtain a 
complementary picture if the continuation is performed while keeping constant another 
characteristic of the breather. In this is done with constant action, here, using 
a similar technique, we have chosen to keep the energy constant, which is physically 
meaningful. 

By fixing the energy we restrict ourselves to the situation which is realized in 
single-molecule experiments. More biologically meaningful would be to keep constant 
not the energy but the temperature. This is however beyond the scope of our 
paper because to carry out this type of investigation it is necessary to introduce into 
equations of motion stochastic forces and solve corresponding Langevin equations. 

The dependence of the breather profile on the parameter fi at constant energy is 
plotted in figure [?]. It can be seen that allowing the breather the freedom of choosing 
its frequency the continuation is possible for much larger values of /i, or, in physical 
terms, there exist localized breathers although not at any frequency. Note that in this 
situation the phonon band is very narrow, and what we have is a nonlinear version of 
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Figure 7. Dependence of the breather profile with hard on-site potential for a 
curved chain with curvature k = 2 on the dipole coupling parameter fi for constant 
energy, E = 0.3855. The frequency changes from 1.2 to 1.46. £ = 0. 



the linear mode described previously, whose frequency is close to the breather's one 
although not identical. 

Although the cases described above are interesting to understand the phenomena, 
the situation with more physical interest is the process of curving the chain. Suppose 
that £ ^ and \i ^ and the breather has its frequency above the phonon band. When 
the chain is curved, the dependency of the breather with respect to the curvature is 
similar as with respect to fi. If we choose to maintain constant the frequency the plot 
of the breather evolution is very similar to figure except for the fact that at the zero 
curvature the breather would have a zig-zag profile which progressively changes to the 
bell profile for large enough curvature. Perhaps more interesting is choice of keeping 
the energy constant, as it is a conserved quantity. The dependence of the breather 
amplitudes and frequency are shown in figure |^. In the figure to the right we can 
see how the linear mode appears and the frequency of the top mode increases until 
almost colliding with the breather frequency. The breather with frequency cjb = 1.2 
no longer exists while there exists a breather which is the nonlinear analogue of the 
top linear mode and its frequency increases as the chain is curved. 

The linear mode is a solution of equation (|J), with W — 0, or, equivalently, with 
almost zero amplitudes for W ^ 0. What we have called its nonlinear analogue, is 
solution of the same equation with \& ^ and has nonzero amplitudes, has a similar 
profile, as can it is shown in figure g[ It has more than a single harmonic in its linear 
spectrum. These extra harmonics are significatively larger for breathers with soft on- 
site potentials, as studied in the next section. In the case studied here, as the linear 
mode is localized, and, therefore, its frequency is isolated from the phonon band, the 
nonlinear analogue can be obtained by substituting ty — > sty in equation (^), and 
changing continuously s from zero to 1 at constant energy. 
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Figure 8. Left: Dependence of the breather profile with hard on-site potential 
for a curved chain on the curvature at constant energy. Right: Plot of the 
breather frequency (dots) and the linear spectrum with respect to the curvature. 
E = 0.3064, e = 0.05, (i = 0.5. 

The observation of figure || shows that for frequencies above the top of the phonon 
band, i.e., slightly below 1.1 for the parameters in the figure, it is only possible to 
continue the breather with the curvature at constant frequency until it approaches 
the top mode frequency curve. Thereafter, it is possible to continue it at constant 
energy, and its frequency follows slightly above the top mode. 

If the chain is curved around the neighbouring site where the initial excitation 
is located, it will switch to the bending point, for even small curvature, when the 
frequency of the top linear mode approaches to its frequency. This happens at 
k = 0.7 for the values corresponding to figure |[ Similar behaviour has been found in 
reference [pi. 

Let us mention that all the breathers described here are linearly stable, i.e. all 
their Floquet eigenvalues have modulus 1. The hard breathers in this system are 
always stable. The properties of the soft ones are very different. 

5. Soft breathers 

We have included the previous section for completeness, but, actually, in biological 
molecules, chemical bonds are thought to be better described by soft potentials. The 
most commonly used are the quartic soft, V(u n ) — hoJ^u^ — the cubic potential, 
V(u n ) = \oJqu\ — \u^, the Lennard- Jones potential, V(u„) = a(l/r 6 — 1/r 12 ), and 
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the Morse potential V(u n ) = Z?(exp(— b u n ) — l) 2 . The last two are known to have 
several nice characteristics: a) they are asymmetric, with a hard part, that describes 
the strong repulsion when two atoms or molecules approach, and a soft part that 
becomes flat, reflecting the weakening of the bond when the molecules get separated 
and eventually unbonded; b) they have moving breathers with stacking coupling 
potentials; c) compared with another frequently used soft potential, the cubic one, 
it does not have an infinite well, which has no physical interpretation and can produce 
anomalous results when simulations are done. We have chosen for presenting our 
results the Morse potential because it is mathematically simpler, but the Lennard- 
Jones potential provides similar results. 

For the straight chain and the normalized Morse potential D = 1/2 and 6 = 1, 
which gives the same frequency ujq = 1 for small oscillations, and a representative 
frequency of the nonlinear excitations lu^ — 0.8 , the breathers are bell-shaped, their 
amplitudes grow with the coupling parameter and they become unstable for e = 0.12. 
The nonlinear mode that produces the instability is an spatially asymmetric one. This 
will be of special interest in the section of moving breathers. 

For stationary breathers, we are interested in values of the parameters where the 
dipole interaction is significant and the nonlinearity is weak. In this situation, typical 
parameters can be e — 0.05, /i = 0.5. For these parameters the frequencies of the 
phonon band, i.e., the continuous part of the linear spectrum, spread from 1 to 1.1. 
A rough measure of the nonlinearity of the excitation is the distance of its frequency 
from the phonon band, which suggest a frequency around 0.8 or 0.9. 

If we start with the straight chain and a breather centered at the bending point 
and we increase the curvature at constant frequency, the same difficulties as with the 
hard breathers arise: the amplitudes of the breathers tend to zero when its frequency 
approaches the bottom linear mode and the continuation becomes impossible. We 
must emphasize this point: in the breather literature path continuation is usually 
done at constant frequency, which is a natural consequence of the theorems of existence 
JTl| from the anticontinuous limit. At the anticontinuous limit there exists freedom 
to choose the frequency of the isolated oscillators and therefore at low coupling too. 
This is not the case if the system is in presence of strongly localized modes: there are 
no localized excitations at every frequency. This characteristic opens the possibility 
of spectroscopic analysis to check the existence of breathers in biomolecules and 
eventually to measure curvatures. 

From the mathematical point of view, if the frequency is not fixed, we need to fix 
another quantity, and the chosen one, as previously, is the energy. Thus, the variation 
of the breather with the curvature mimics an adiabatic process of bending. Figure 
left, shows the dependence of the breather profile with u>b = 0.8 at k = on the 
curvature, changing its shape from a bell profile to a zig-zag one, until its becomes 
almost invariable. While this process takes place, its frequency changes, which is 
shown in figure ^-right. The breather becomes unstable at k w 0.8. The frequencies 
of another breather with lower energy and higher frequency is also represented in 
this figure. The higher the frequency the larger becomes the curvature, for which 
instability is produced. 

The exact stability analysis has been done calculating the Floquet eigenvalues, 
but we prefer to plot here the breather frequency, the linear spectrum and 1.5 times 
the breather frequency. This gives a much more intuitive physical picture. 

The reason for plotting this frequency is the following. When the curvature is 
increased, the frequencies {tUi\ of the internal modes, i.e., the small perturbation of 
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Figure 9. Left: Dependence of the breather profile with soft, on-site potential for 
a curved chain on the curvature at constant energy E=0.35. Right: Plot of the 
breather frequencies for two energies, E=0.35 (big dots) and E=0.20 (dashed), 
and 1.5 times the breather frequencies (small dots) and the linear spectrum 
(continuous lines) with respect to the curvature. Parameters: e = 0.05, fi = 0.5. 

the breather, are very similar to the ones shown in figure |], as commented in section [| 
The biggest difference is that the bottom mode (with soft on-site potential, and the 
top one for hard on-site potential) is not there, because this mode is substituted by the 
breather itself. The corresponding Floquet multipliers are {A; = exp(± 27rw;/wb)}, 
{9i = ±2nuji/uJb} called the Floquet arguments. For a instability to be produced, 
a collision of two multipliers has two take place with the condition that they have 
the same Krein signature (See reference [|l3| for details on this subject). The positive 
(negative) Floquet arguments have all of them the same positive (negative) signature 
and cannot therefore bring about instabilities until they collide at (1 + Oi), named 
harmonic instabilities, or (—1+0 i), subharmonic instabilities, or between them, named 
oscillatory instabilities or Krein crunches. The first one happens when 9i = 2 7rn, 
n G N, for some i, that is, u>i = nuj^, or in other words a linear mode resonates with 
the breather frequency or its harmonics. Figure |^ shows that this kind of bifurcation 
does not happen here. 

The subharmonic instabilities can occur when 0i = n + 2n7r,, for some i, i.e., 
oji = Ldh/2 + nuib- Within the range of parameters studied, this can only occur when 
some u>i equals 1.5 ujb- Therefore, subharmonic instabilities appear in our system, when 



Breathers in curved alpha-helix 



13 



some frequency of the internal modes (very close to the linear ones, as commented 
above), collide with the 1.5 Wb- This is easily seen plotting simultaneously the 1.5wb 
curve and the linear frequencies. 

Figure |^— right shows the reason of this instability: the curve 1.5 Wb intersects 
with the frequency curve of the top linear mode, bringing about the instability. 
Nevertheless, the latter is a localized mode, centered at the bending point and the 
localization persists although the excitation has now superimposed a perturbation 
with twice the breather period. Within the range of parameters shown there are also 
a Krein crunches, related with perturbations of the frequency, but also preserving 
the localization. This happens after the Floquet multipliers of the phonon band have 
crossed at (— 1 + Oi), or, in other words, their frequency have become larger than 
1.5 u>b- The eigenvalues of the internal linear modes cross with other and do small 
excursions outside the unit circle. This is also evident, as the curve 1.5 Wb enters the 
phonon band. The same figure shows the frequency of a breather with higher energy 
corresponding to a frequency Wb = 0.9 at k = 0. Note that both breathers tend to 
the same frequency as the curvature increases and that breathers with energy slightly 
below E = 0.35 avoid the secondary bifurcations, because the 1.5 Wb curve does not 
intersect the phonon band. 

The apparent intersection of the breather frequency with the bottom linear mode 
is misleading because the latter does not exist as a perturbation of the breather, being 
this one its nonlinear analogue. The breather is not, however, linear, as its first AC 
harmonic is of the order of the DC one, and the second, about one third of it. 



6. Moving breathers 

In this section we explain briefly the behaviour of a moving breather in a curved region 
of the alpha-helix. 

In order to move a breather, we apply the marginal mode method [ [17| . It basically 
consists in adding a spatially asymmetric mode (the pinning mode) to the breather 
velocity. As is shown in p0[ , the stacking coupling of the breather must be strong 
enough in order that it can be moved. 

The scenario is similar to the observed in DNA chains Jl(| and biopolymers 
i.e., for a fixed value of the curvature, there exists a critical value of the velocity below 
which the breather is reflected when it reaches the bending point and, above which, the 
breather crosses through it. Analogously, for a fixed value of the velocity, the breather 
crosses the bending point as long as the curvature is smaller than a critical value. As 
the curvature increases, the breather spends more and more time at the bending point 
and eventually is reflected. These two behaviours are shown in figures |l^ and [ll[ In 
other words, the moving breather in a curved alpha-helix chain behaves as a particle 
in a potential barrier. 

In reference pl| an hypothesis is introduced about the existence of trapping in 
inhomogeneous Klein-Gordon lattices. According to it, trapping of breathers does not 
occur when there exists a linear localized mode with a profile (vibration pattern or 
wave vector) different from the stationary breather one. When this occurs, the moving 
breather behaves as a particle in a potential barrier and is never trapped. As Figure [| 
shows, there is a localized mode (the top mode) of opposite profile to the breather, and 
the case exposed in the reference cited above is the same as in the curved alpha-helix. 
The simulations confirm that the trapping hypothesis pl[ | holds in our system, and 
moving breathers cannot be trapped. 
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We still have not been able to find a counterexample, and this hypothesis seems 
generic for Klein-Gordon systems. However, it has no mathematical proof and should 
be checked as much as possible to achieve, at least, an inductive confirmation. 

7. Conclusions 

We have studied the existence and properties of stationary and moving breathers in 
a model with dipole moments parallel to a curved chain of oscillators. Its interest is 
twofold. On the one hand it is a model for a channel of amino acids along an alpha- 
helix protein, where the internal degrees of freedom are the Amida-I vibrations. In 
this aspect the objective of our work is to explore the role of the bendings as places 
where the energy is stored and to investigate whether moving excitations can travel 
across a bending point or not. On the other hand the theoretical interest comes from 
it being a model that completes previous works on bending chains. The specifics of 
the model are: a) the dipoles are oriented along the chain instead of perpendicular to 
it; b) it is sufficient to take in account nearest neighbours in order that the system 
can feel the shape of the chain. 

In previous models [?], |l(J the shape is felt due to the change between 
the distances of the oscillators, and, therefore, they need long-range interaction in 
inextensiblc chains. In our model the shape is felt due to the change of the angles 
between dipole moments. In this way, we expect to help to complete the study of 
systems in curved chains and approach to a generalized description. 

The most important general fact, overlooked in the breather literature to our 
knowledge, is that the presence of strongly localized models around the bending point 
has such decisive a role: the breathers cannot be chosen at any frequency when the 
curvature is increased. Their interaction with the linear localized modes brings about 
frequencies approximately determined, and dependent on the curvature. 

Consequently, we have developed a variant of the existing techniques, described 



in Appendix B to obtain breathers at constant energy. Therefore, the continuation 
of the breathers when the curvature is increased, models an adiabatic process of 
bending. More biologically meaningful would have been to model this process at 
constant temperature, but this is outside the scope of the techniques used in this 
paper, which are related more to single molecule experiments. 

The adiabatic process of curving the chain transforms the initial breather into a 
nonlinear analogue to the top or bottom linear mode according to the on-site potential. 
The corresponding linear mode, therefore, do not exist in presence of the breather 
and cannot produce instabilities. The consequence is that hard breathers in the bent 
chain are stable and soft breathers experience subharmonic instabilities, that conserve, 
however, the localization. 

We also have shown that up to curvatures relatively high, moving breathers can 
travel across the bending point spending some time at it. They cannot, however be, 
trapped, which confirms recent developments in the field j2j| . 

The physical suggestion is that breathers may be a means of storage and transport 
of energy in proteins, due to and in spite of their complex structure in the space. 
Certainly, our model is a rough description of such complex a system, yet the main 
conclusion of the present and previous studies could be that the properties of breathers 
in different curved chains are quite common. We also think that they might provide a 
means to detect experimentally breathers at the bending points, showing absorbtion 
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Figure 10. A breather travelling across the bending point. A small increase of 
the curvature will produce a reflection 



of characteristic frequencies and their harmonics, which depend on the curvature. We 
are, presently, in contact with experimental groups to explore its feasibility. 

Appendix A. Details of the model 

As said in Section [l], we describe the position of each amino acid by a vector r„, n 
being an index. The distance between amino acids is supposed to be a constant a, 
i.e. , a = |r„ + i — r„|. The direction of the dipole moments are given by unit vectors 
t„ = (r, i+ i — r„)/|r, i+ i — r„| and the moments themselves are given by p n = p n t„. 
The interaction energy between two neighbouring dipoles n and n + 1 is given by 

Pn+lPn Q Pn(r„+l ~ r n) ' Pn+l(l*n+l ~ r n) 



U„ n +l — 



|r n+ i-r„| J |r„ + i-r„| D 

tntn+i „t n (r„ +1 — r„) • t n+1 (r„ +1 — r„) 



Pn+lPn 



|r, i+ i-r„| 3 |r„ + i-r„| 5 / (A.l) 



Pn+lPn o 3tt 1 

a J a 



= (Pn+iPn/a 3 ) (-2t n+ it„) = (p n+ ip n /a 3 ) ((t n+ i - t„) 2 - 2) 

Suppose that the equilibrium value of p n = po, Vn, then we can express its 
non equilibrium value as p n = po + qu n , where u n represent the stretching from the 



equilibrium positions. Substitution into Equation A.l leads to 



2 2 

U n ,„+l = f/°,„+i+^, n+ i+( -2^ + ^(tn+l - t„) 2 ) U n + 1 U„ (A.2) 
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Figure 11. A breather being reflected at the bending point. 



The first two terms in this equation stand for constant and linear terms, therefore, 
they will not appear in the Hamiltonian, as its first derivatives with respect to the 
variables u n at equilibrium, i.e., u n = must be zero. We represent the remaining 
term by H n>n+1 . 

Suppose that the Hamiltonian of the system is given by 

H = J2 \ul + \a 2 u 2 n + + ^K(u n+1 - u n f + H n , n+l , (A.3) 

n 

where a and K are positive constants and f(M„) is the nonlinear part of the on-site 
potential. Rearranging the terms we obtain 



*(«n) 



( K + 22 3)( U «+1 _ U «) 2 + l(*I>+l - t n ) 2 U n+ lU r , 



2 V a 

Thus, the Hamiltonian can be written as 



(A.4) 



H = ^ + ^ w u n + *(«n) + ^K+l - U n ) 2 + (J,(t n +1 - t n ) 2 U n+1 U n , 

n 

with u>o = 1, e = (K + 2 q 2 /a 3 )/a' 2 , fj, = q 2 /(a 3 a' 2 ), where a' 2 = a 2 - 4q 2 /a 3 , with 
the time rescaled t — * a'i and the redefinition of &(u n ) — > 4 , (M n )/a' 2 . The initial 
Hamiltonian has been also divided by a' 2 , 

We keep the term wq in the formulae in spite of being 1 as a reference because 
its meaning is the frequency at very low amplitude when the oscillators are isolated. 
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The corresponding dynamical equations are 



U„ + LO^Un + *'(u n ) + e(2 u n - u n+ i - u n ) 



+ /u((t„+i - t„) 2 u„+i + (t„_i - t„) 2 M«-i) = 0, (A.5) 
which are written in a simplified notation in Equation |^. 

Appendix B. Breathers with constant energy 

To obtain breathers with constant energy, we use a variant of the method in the time- 
Fourier space Jl^, [23|l . First, we obtain breathers by path continuation at constant 
frequency from the anticontinuous limit using the Newton method as described in the 
references above. Thereafter, the path continuation is performed at constant energy. 
In this case our variant applies and we give below the details needed to obtain the 
corresponding Jacobian used by the Newton method. It is also a variant of the method 
developed in 22 at constant action. 



The Hamiltonian of the system can be written: 

H=±Y. il l + W(u,n), (B.l) 

n 

where u = (u±, . . . , ujy), and W{u, n) represents the potential energy, sum of the 
on-site and coupling energies; n represents any parameter, as the curvature in this 

dW 

article. The dynamical equations are given by u n + — — = 0. Time-reversible periodic 

ou n 

solutions with frequency Uh are given by a truncated Fourier series: 

k m 

u n(t) = J^(2 - 4,o) z k , n cos(kuj h t) (B.2) 

fc=0 

The use of functions with different frequencies is inconvenient and can be avoided by 
the change of variable t — u^t, which leads to the Hamiltonian: 

H=\cjIY / <" + W(u,k), (B.3) 

n 

where / represent the derivative with respect to t. The corresponding dynamical 
equations are: 

f n (u,cj h ,n) =wg< + f^ =0 (B.4) 

In this equations LUh enters as a parameter and the functions «„, with frequency unity 
in the new time variable, are given by: 

u n{t) = ^(2 - 5 k ,o)z k ,n cos(fcf) (B.5) 

k=0 

Therefore, the functions /„ can be seen as 27r-periodic functions of t due to their 
dependence on u(t) and u"it). 

The c omp onents of the cosine discrete Fourier transform of the dynamical 
equations (|B.4| ) are given by: 

Fk[fn] = -k 2 culz k , n + F k [^(u(t))} , (B.6) 
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where F k [g(t)] represent the k-th term in the cosine discrete Fourier transform of a 2n 
periodic function g(t), given by: 

k 

x-, r 1 /„ r \ / 2wm , , 2Trkm . . . 

Fk[9(t)} = — £ (2 - J TO)0 ) g( — ; J cos( o; _ ; J (B.7) 

m— 



l 



2 /c^j 



In this way, finding a breather is reduced to solve the (k m + 1) x N equations (B.6). 
As a breather solution is determined by the same number of Fourier coefficients z k ^ n 
and the frequency LOh, we need and extra equation. If we are interested in obtaining 
solutions with given energy E, this is given by H — E = 0, which for time-reversible 
solutions (u' n (0) — 0) reduces to 

F H = W(u(0),k)-E= 

We denote by z and F the column matrices ({#fc,7i}, Wb) and ({Fk [fn]}> Fh)> 
respectively, with the index fc running faster. In order to use the Newton method for 
path continuation with the parameter k we need the Jacobian: 

dF k [f n 

[dF k [f r r 



dF 

dz 



I 9zk', n 
dF H 

dZk'n 



duib 

dFn 
duj h 



(B. 



The left upper block is the Jacobian at constant frequency, its elements given by: 



dF k [f n 

dZk'n' 



k 2 UjlS k ,k'S n ,n' + 



2 - 5 ,, 



Fi 



d 2 W 



k+k> 



du n du r , 



+ F\ k _ k ,\ 



d 2 W 

dUndUn 



(B.9) 



The elements in the last column are given by: 



dF k [f n 



= -k 2 U>bZk,: 



dF, 



H 



dW(u(0)) 



The remaining elements in the last row are: 
dF H _ dW({u n (0)}) _ aW{{Y^o(2-6i, )zi, m }) 



duj h 



(B.10) 



dz k > 



dz k > 



dz k , 



dW \ 

du n , ) t=Q 



(B.ll) 



2=0 
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